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The electronic quasiparticle spectrum of a ferromagnetic film is 
investigated within the framework of the s-f model. Starting from 
the exact solvable case of a single electron in an otherwise empty 
conduction band being exchange coupled to a ferromagnetically sat- 
urated localized spin system we extend the theory to finite tempera- 
tures. Our approach is a moment-conserving decoupling procedure 
for suitable defined Green functions. The theory for finite tempera- 
tures evolves continuously from the exact limiting case. The restric- 
tion to zero conduction band occupation may be regarded as a proper 
model description for ferromagnetic semiconductors like EuO and 
EuS. Evaluating the theory for a simple cubic film cut parallel to the 
(100) crystal plane, we find some marked correlation effects which 
depend on the spin of the test electron, on the exchange coupling, 
and on the temperature of the local-moment system. 



I. INTRODUCTION 

Since the mid-70's there hasJheen a growing theoretical in- 
terest in the s-f (or s-d) modelErcl. The model works for ma- 
terials which exhibit local-moment magnetism: for magnetic 
semiconductors like the europium chalcogenides EuX (X=0, 
S, Se, Te)H and for metallic local-moment systems such as Gd, 
Tb, and Dyu. In the local-moment magnets, the electronic 
and the magnetic properties are caused by different groups 
of electrons. Whereas the electronic properties like electrical 
conductivity are borne by itinerant electrons in rather broad 
bands, e.g. 6s, 5d for Gd, the magnetism is due to a strongly 
localized partially filled 4/-shell. In the case of Gd and Eu 
compounds the 4/-shell is exactly half-filled and, because of 
Hund's rules, has its maximal magnetic moment of S = |. 

Many characteristics of the local-moment systems may be 
explained by a correlation between the localized magnetic 
states and the itinerant electrons. In the s-f model this correla- 
tion is represented by an intra-atomic exchange interaction. 
The difference between the s-f model and the well-known 
Kondo lattice modela is that in the former the exchange inter- 
action is ferromagnetic, favouring parallel alignment of itin- 
erant electrons and local-moments, whereas in latter it is an- 
tiferromagnetic. That is why the s-f model has been recently 
more and mote, often referred to as the Ferromagnetic Kondo 
lattice modeWa. 

The second aspect of this paper is that of reduced dimen- 
sionality. Magnetic phenomena at surfaces and in thin films 
attract broad attention both theoretically and experimentally 
due to the question of phase transitions and the variation of 
magnetk and electronic properties in dimensionally reduced 
systemsETtil. One of the most remarkable examples of the 



outstanding magnetic properties at surfaces is the existence 
of magnetically ordered surfaces at temperatures where the 
bulk material is paramagnetic. This effect has been, first 
documented for Gd(0001) surfaces by Weller et al.tB and 
since then been measured by difiejrenJ: groups using a wide 
range of experimental techniquesll3~H In these experiments 
for the difference between the Curie temperature at the sur- 
face, Tc (surface), and the Curie temperature H&f bulk Gd, 
Tc (bulk), values between 17M13 and some 60KLL3 have been 
reported. Succeeding the results for Gd, Tb also was found 
to haye_a higher surface Curie temperature, relative to the 

buMM n 

Contrary to the groups cited above, Donath et altll, using 
spin-resolved photoemission did not find any indication for 
an enhanced surface Curie temperature of Gd(0001) surfaces. 
Other controversially discussed surface properties of Gd in- 
clude theieirmerature dependent behaviour of a Gd(0001) sur- 
face statecirtJ which is supposed to play an important role in 
the interplay between electronic structure and magnetism. A 
thorough account on the surface magnetism of the lanthanides 
has been recently given by Dowben et alt3. 

It is not only the dimensionally reduced Gd which is of in- 
terest here, even bulk Gd is far from being-.completely un- 
derstood. In an earlier study Nolting et alx3 have predicted 
that the a priori non-magnetic (5£/,6s)-conduction and va- 
lence bands should exhibit a marked non-uniform magnetic 
response at different positions in the Brillouin zone and for 
different subbands.. Weakly correlated (s-like) dispersions 
show a Stoner-like T-dependence of the exchange splitting. 
On the other hand, stronger correlated (cZ-like) dispersions 
split below Tq into four branches, two for each spin direction. 
Their T-dependence mainly concerns the spectral weights of 
the quasiparticle peaks and not so much the energy positions. 
Consequently, an exchange caused splitting remains even for 
T > T c . This may be the reason for the fact that the ex- 
perimental situation is controversial. Kim et al£B found a 
T-dependent spin splitting of occupied conduction electron 
states, which collapses in a Stoner-like fashion for T — ► Tc- 
From photoemission experiments, Li et a/.c3 E3 conclude that 
the exchange splitting must be wave-vector dependent, col- 
lapsing for some k values, while for other no collapse occurs 
as a function of increasing temperature. This fairly compli- 
cated temperature behaviour in the bulk-material must be ex- 
pected for Gd-films, too. 

It is not at all a trivial task to perform an electronic structure 
calculation for a ferromagnetic local-moment film in such a 
manner as to realistically incorporate correlation effects. In 
a previous papetca we proposed a simplified model which 
allowed us to exactly calculate the electronic structure of a 
model film in the limiting case of ferromagnetic saturation and 
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empty conduction band, n = 0. This case is applicable to a 
film of a ferromagnetic semiconductor such as EuO, EuS at 
T = 0. Its significance arises from the fact that all relevant 
correlation effects which are found or expected tojaccuj^at fi- 
nite band occupations and arbitrary temperaturescSEMJ^Ldo 
already appear in this rigorously tractable special caseErtS. 

In this paper we extend the T = special case to finite tem- 
peratures. As a result we will calculate the electronic structure 
of a local-moment film with a single-electron in an otherwise 
empty conduction band within the whole temperature range 
from T = and T = T c . 

In the next section we present the model and define the 
corresponding many-body problem. S ubse quently the model 
is evaluated in two steps, in s ection II A for the electronic 
subsystem, and in section [IB for the local-moment system. 
Section 



III 



is devoted to a detailed discussion of the results 
obtained for different film thicknesses and various exchange 
couplings and temperatures. Comprehensive conclusions with 
an outlook on the possible application of the model to real 
substances and on the evaluation of temperature dependent 
surface states in section [iv| complete the paper. 



II. THEORETICAL MODEL 

We investigate a film consisting of n equivalent layers par- 
allel to the surface of the film. Each lattice site of the film 
is indicated by a greek letter a, ft, 7, ... denoting the layer 
index and a latin letter i, j, k, . . . numbering the sites within 
a given layer. Each layer possesses two-dimensional trans- 
lational symmetry. Accordingly, the thermodynamic average 
of any site dependent operator Ai a depends only on the layer 
index a: 



The complete s-f model Hamiltonian 

TL = TL S + Hf + Tt s f 
consists of three parts. The first 

ijafi 



(1) 



(2) 



(3) 



describes the itinerant conduction electrons as s-electrons. 
cf aa and Cj/3 a are, respectively the creation and annihilation 
operators of an electron with the spin a at the lattice site Ri a . 
T^j are the hopping integrals. 

Each lattice site Ri a is occupied by a localized magnetic 
moment, represented by a spin operator SjQ,. These localized 
moments are exchange coupled expressed by the Heisenberg 
Hamiltonian: 



Mermin-Wagner theoremcS there is no solution showing col- 
lective magnetic order at finite temperature, T > 0. To avoid 
this obstacle we have chosen an extended Hamiltonian for the 
localized moments, 



rtt = rif + tiA = rlf 



(5) 



which additionally to the Heisenberg Hamiltonian, Tit fea- 
tures a single-ion anisotropy term Ha- Dq is the according 
anisotropy constant, which is typically smaller by some or- 
ders of magnitude than the Heisenberg exchange interaction, 

Do « Jf. 

The distinguishing feature of the s-f model is an intra- 
atomic exchange between the conduction electrons and the 
localized/-spins, 



H s f = 



(6) 



Here, J is the s-/exchange interaction and a ia is the Pauli spin 
operator of the conduction band electrons. For the materials 
we are interested in the s-f coupling is positive (J > 0). In 
the case where J < the model Hamiltonian (^) is that of the 
so-called Kondo lattice. Using the second-quantized form of 
ai a and the abbreviations 



q 



z Ul) 



±1, 



the s-f Hamiltonian can be written as 



H 



-E 



(7) 



(8) 



The most decisive part of the s-f Hamiltonian (||) is the second 
term, which describes spin exchange processes between the 
conduction electrons (fy and the localized moments (Q). 

In general, the alignment of the localized moments will be 
influenced by the s-f interaction, which can mediate an indi- 
rect interaction (RKKY) via the occupied conduction bandEJ. 
However, here we are interested in the electronic quasiparti- 
cle spectrum of a ferromagnetic semiconductor according to 
a simple test electron in an otherwise empty conduction band. 
In this case, the localized spin state cannot be affected by the 
s-f interaction. Furthermore, one knows from experiment that 
typical Heisenberg exchange integrals are smaller by some or- 
ders of magnitudes than their s-f counterparts. In this respect, 
it seems appropriate to neglect the Heisenberg exchange in- 
tegrals for the calculation of the electronic properties of the 
system. Accordingly, the Hamiltonian ^) can be split into an 
electronic, H s + W s /, and a magnetic part, 7i*p which can be 
solved separately. 



ija/3 



(4) 



where are the exchange integrals. The problem with the 
simple Heisenberg model in the form (0) is that due to the 



A. The electronic subsystem 

Starting from the Hamiltonian of the electronic subsystem, 
H* = H S +H sf , (9) 
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all physical relevant information of the system can be derived 
from the retarded single-electron Green function: 



-l j dl v 



-±Et 



(10) 



Here and in what follows [.,.]+ ([.,.]_) is the anticommuta- 
tor (commutator). Conform to the two-dimensional transla- 
tional symmetry, we perform a Fourier transformation within 
the layers of the film, 

G ^{E) = l^e ik ( R -^)G«e(i?), (11) 



where N is the number of sites per layer, k is an in-plane 
wavevector from the first 2D-Brillouin zone of the layers 
and Ri represents the in-plane part of the position vector, 
Ri a = Ri + r a . From Eq. (|ll]) we get the local spectral den- 
sity by 



(12) 



which is directly related to observable quantities within an- 
gle and spin resolved direct and inverse photoemission exper- 
iments. Finally, the wave-vector summation of S^(E) yields 
the layer-dependent (local) quasiparticle density of states: 

k 

In the following discussion all results will be interpreted in 
terms of the spectral density ( [12] ) and the local density of states 

®" ^ 

For the solution of the many -body problem posed by Eq. (|9J) 

we write down the equation of motion of the single-electron 
Green function ( |To[ ) 



in ji 



(14) 



where 6"? = 5 a p$ij- The formal solution of Eq. ( Jl4| ) can be 
found by introducing the self-energy M^(E), 



= ^M1(E)G^(£), (15) 



which contains all information about the correlations between 
the conduction band and localized moments. After combining 
Eqs. ( p"4| ) and ( p"5| ) and performing a two-dimensional Fourier 
transform we see that the formal solution of Eq. (|l4j) is given 
by 



Gker(-E') = H(E1— Tk — Mk<j (E)) 



(16) 



where I represents the (n x n) identity matrix and where the 
matrices G] i(7 (E), Tk, and Mko-(-E) have as elements the 
layer-dependent functions G^(E), T^, and M^(E), re- 
spectively. 

To explicitly get the self-energy in Eq. ( p"5[ ) we evaluate the 
Green function 



— ( ' ■? r aQ/3 

" 2 V Hj° 



tpaa/3 
iija 



Here the two higher Green functions, 



= ((^ fc7 -,;c+ CT )) 



(17) 

(18) 
(19) 



originate form the two terms of the s-f Hamiltonian (|8|) 
and will be referred to as the Ising and the Spin-flip 
function, respectively. Considering the equations of mo- 
tion for these two Green functions we encounter the 
two higher Green functions ({[Sf a c kl<T , T~L s f]-\ c%J)e and 

({[^ia Ck-y-v, ~H s f]-; ctpj) £• Since we consider an empty 
conduction band the thermodynamic average in the Green 
functions has to be computed with the electron vacuum state 
\n = 0). From the definition of the s-f Hamiltonian (|J) we 
then see that (n = 0| H s f = and, accordingly, 



[S- a ,Hsf}_c kl< T;> 



0. 



s/J _ ca_, 7 _ ct ; Cj j3c 



0. 



Hence, for the equations of motion of the Ising and the Spin- 
flip function we get 

= n (s z a ) qf + ({sf a [c k ^,n sf }_ ; c + a )) E (20) 



nifi. 



(21) 



On the right-hand side of these equations appear further higher 
Green functions which prevent a direct solution and require 
an approximative treatment. The treatment is different for 
the non-diagonal terms, (i,a) ^ (fc,7) and for the diago- 
nal terms, (i,a) = (k, 7). In the first case we use a self- 
consistent so-called self-energy approach which results in a 
decoupling of the equations of motion. For the diagonal terms, 
(i, a) — (k, 7), this approach is replaced by a moment tech- 
nique which takes the local correlations better into account. 

a. Non-diagonal terms (i,ot) ^ (k,j): The definition of 
the self-energy ( p"5| ) formally corresponds to the substitution 



?mtr 



(E)c n 



(22) 
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within the brackets of the Green function. The inspection of 
the spectral decomposition of the two functions in Eq. ( |l5| ) 
reveals that both, (({c iaa ,H s f}-; c^}) E and ((<w; cf ffa .)) E , 
have the same pole structure and can differ only by the spectral 
weights of their poles. The equality of both sides in Eq. ( |l5| ) is 
installed by the self-energy components M?P(E). Inspecting 
now the spectral representations of the two Green functions 



({S ia a [ck^,n s f}-;cj jE and ((S la a c kj<T ; c+ J E we no- 
tice that the additional spin operator Sr CT selects for both only 
those poles of the original Green functions without spin op- 
erator which are connected with a spin-flip of the electron. 
Hence, the poles of these two functions build subset of the 
poles of the two Green functions from Eq. (^5]) and are, there- 
fore, identical to each other. Again, only the weights of the 
poles can differ. In analogy to Eqs. ( |15| ) and (|2|) we now 
propose to use the plausible ansatz 



S ia i C kl-a,'H s f}_;- 



(23) 



A similar reasoning can be used for: 



E 

m ft 



(24) 



with the difference that here the additional spin operator Sf a 
does not change the original pole structure but merely modi- 
fies the spectral weig hts of the poles. On the right-hand sides 
of Eqs. ( P3[ ) and (£4j) we find the already known Spin-flip and 
Ising function, respectively. Hence, for (i,a) ^ {k,j), the 
Eqs. @, (|n|), ©, (0), (||), and (|f) build a closed sys- 
tem . 

b. Diagonal elements (i, a) — (k,j): We start with the 
explicit evaluation of the higher Green functions on the right- 
hand sides of Eqs. ( pp| ) and (pT|). For Eq. ( j2l| ) we get, for 

(i,a) = (k, 7), 



As Eq. (£5|), the above relation is still exact. To get a close 
system of equations we are left with the determination of the 



functions F" t "f(E) and F£?£(E). Both fulfil exact relations 
which will be used to derive satisfying approximations. For 
spin S — h we find for all temperatures: 



acxfS 1 



paa/3 



\z a hP^{E), 
±h 2 G»f a (E)- Z(T hr%°f 



(E). 



(28a) 
(28b) 



On the other hand, in the case of ferromagnetic saturation, 
{S a ) = S, it holds for arbitrary spin: 



&(E) 



T=0 



K{{S-\) + \z a )F%*[E), 



F^{E) 



T=0 



= H 2 SGf a (E) 



The exact limiting cases ( |28[ ) and 
structures: 



(29a) 
(29b) 
suggest the general 



it j a 

z a nr^{E) 



(30a) 
(30b) 



For the five Green functions of the type ({A\ B)) b in Eqs. (^ 
we can calculate the spectral moments, 



1V1 AB — 



ifi|) [.!(/). #(())]_ 



t=o 



where n = 1,2,.... Because of the equivalent relation 



M (n) - - 

1U AB — 



dE E n lm ((A; B)) 



E • 



(31) 



(32) 



the moments can be used to fix the coefficients and AS/ 



in Eqs. 
we get 



After tedious but straightforward calculations, 



S ia [Cia-cr, 7~L s f]_ ; Cjp^ 



where we have abbreviated 



Fuji (F) = \\S ia a S° a c iaa -, c^pcr 



(25) 

(26a) 
(26b) 



The analogous evaluation of the higher Green function in Eq. 
( pc| ) does not require any further higher Green functions, be- 
cause it can be expressed in terms of already known Green 
functions: 

{{Si a \ c ^a , H sf ] _ ; cfpj) E + z a {{S^J [c iQ _ CT , H sf ]- ; c+ ja )) E 

= ™{F^{E) + z a E^(E) - z a hS(S + l)Gf a (E)). 

(27) 



f^aa — 0, Kq(| — {S a a S^) ^ai (^a)i 



_ {S a a SgSq 



Z<r{S a a Sa) 



(33) 



, (2) _ (Sa^SgSa} - (S^iSa^Sg) 

<(^) 2 > - {S^ 

The coefficients are determined by /-spin correlation func- 
tions, which will be determined at a later stage. 

The Eqs. @, @, @, @, ©, ©, and 

( |33| ) represent a closed system, which can be solved self- 
consistently. Before proceeding we assume that the self- 
energy from Eq. (|l5|) is a local entity 



Mf a {E) = S aP M«(E). 



(34) 



The reason for the k independence can be traced back to the 
neglect of magnon energiesta. The restriction to the diago- 
nal elements in the greek indices denoting the layers is in that 
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sense the transfer of the k independence in the case of three 
dimensionso to the film geometries discussed in this paper. 
Furthermore one can show that the assumption ( [34] ) is not 
necessary for the following calculations but merely drastically 
simplifies them. 

We can now use Eqs. (||)-(||), (§|), (|o|), and S) to eval- 
uate the Ising and the Spin-flip functions in Eqs. (|20|) and (j2l|). 
As the result we get the Fourier transformed Ising and Spin- 
flip functions. According to Eq. (ftjn we can restrict our atten- 
tion to the diagonal elements r°^{E) and F^(E). After 
subsequent q-summation we eventually get, using Eqs. ( |T6| ) 
and@): 



^E r z£ = ^<^>< - gz{m? e r Z 





kqtr 



J [z a {n<Zl-tfS{S + l))VNGt 



(35a) 



a0 
krr 



q 

= - G 0-<r{ M -a £ F kq/ + { ("£ 

q q 

-*a£ + (35b) 

q q 
where we have introduced 

G o"(^) = ^E G k"(^)- (36) 



The set of equations (|35|) can be solved to express the 

SUmS Eq r kqa( E ) aIld Eq F kqa ( E ) in tel ' mS ° f the single- 

electron Green function G^(E). However, by inspecting 
Eqs. ( p5| ) we see that these expressions will still contain the 
layer and spin-dependent self-energy M"(E). 

To solve this problem we combine Eqs. ( p"5j ) and (jl^) and 
get, after Fourier transformation: 



MZGg = - 



J 



2V N 



E^+E^)- (3 7 ) 

q q / 



Combining this equation with the results obtained for 

Eq^kqf ( E ) and Eq^kqf ( E ) from Eqs. © we even- 
tually get an implicit set of equations for the layer and spin- 
dependent electronic self-energy: 



(38) 



where the numerator and the denominator, respectively, are 
given by 

Z a = z„H 2 {S z a ) + J -{(k<£1 - H 2 S(S+1))GX 



-((AW + Ai 2 ) + z a mt a )h{S a ) + Hk^)G^} 

+ ^{z a h 2 S{S + 1)(A« + A£ 2 j + z a m a _ a ) 
+^(m a a - m« a ))GZG a %, (39a) 

K = n2 -^{( h + Z ° X «« + <) G o° + (^ A - 

+m%)Gfc r } + ^{(m« + h)(m% + z CT A«) 
+z„\W(m% + h))GZGZ%. (39b) 



The implicit set of equations (|38j) and (|39j) now enables us 
to self-consistently evaluate the self-energy of the system pro- 
vided that the /-spin correlation functions from Eqs. ([33]) are 
known. These will be evaluated in the next section. 



B. The local-moment system 

The system of localized /-moments is described by the 
extended Heisenberg Hamiltonian (|j) which we write down 
again for convenience: 



H* f = - e j ij s »« s i/3 - A) yz( s i< 

ijctf3 id 



(40) 



Here we want to stress once more that the single-ion 
anisotropy constant Dq is small compared to the Heisenberg 
exchange interaction, Dq <C J%f . By defining the magnon 
Green function 



D°f(E) = ((st a ;Srp 



(41) 



we can calculate the /-spin correlation functions by evaluating 
the equation of motion 



EDtf{E) = 2h 2 5 aP {S z a ) 



,S 



30 



(42) 



The evaluation of this equation of motion involves the de- 
coupling of the higher Green functions on its right-hand side, 
originating from the Heisenberg term. Hf, and the anisotropy 
term, Ha, using the Random Phase Approximation (RPA) and 
a decoupling proposed by LinescJ, respectively_The details of 
the calculation can be found in a previous papeio. For brevity 
we restrict ourselves to present here only the results. For the 
layer-dependent magnetizations of the /-spin system we get 

(S*) (1 + ^ a ) 2S+1 (S - Va ) + vl s +i(S + l + ^ a ) 



fl s+1 - (1 



fa) 



2S+1 



where 



fa = T7 



\ ~" X aa7 (k) 
k 7 



(43) 



(44) 
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where, again, N is the number of sites per layer and (3 = -^p- 
The summation 2~^ 7 in Eq. (44) runs over the n poles B 7 (k) 



of the Green function DZ (E) and the Xaa-y (k) is the weight 
of the 7'th pole in the diagonal element of the Green function 
D^ a (E). The poles and the weights can be calculated from 
the solution of Eq. (42): 



Df(E) = 2h l 



( (Sf) 



V o 







(S Z n ) 



(EI-AY 



(45) 



with 

(A) a/3 

n 



a/3, 



(D <S> a + 2 V J^(S z ))8 a p - 2J^(S z a ). (46) 
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The $ Q come from the decoupling of the higher Green func- 



tion on the right-hand side of Eq. ( |42| ) which originates, 
the anisotropy Hamiltonian TLa according to Lines! 
are given by: 



$a = 



2h 2 S(S + 1) - 3h(S z a )(l + 2^ a ) 
(S z ) 



rom 
and 



(47) 



Having obtained the layer-dependent magnetizations ( |43[ ) 
and the ip a we can now express all the other /-spin correlation 
functions appearing in Eqs. (33 ) via the relations 



((s z a ) 2 ) 



H 2 S(S+l)(S z a )(l + 2^ a ), 
h 3 S(S+l)^ a + H 2 (S z a )(S(S- 
-h{{S z a ) 2 ){\ + 2^ a ), 



and the general spin-operator equality 

SLS£ = t?S{S + 1) + z a hS z a - (S z a ) 2 . 



(48a) 
(48b) 

<Pa) 

(48c) 



(49) 



Mediated by Eqs. (|33[), Q38J), and (|39j), the /-spin correlation 
functions contain the whole temperature dependence of the 
electronic subsystem (ph. 



III. RESULTS 

We have evaluated our theory for a film with simple cubic 
(s.c.) structure consisting of n layers parallel to the (100)- 
plane of the crystal. The electron hopping and the Heisenberg 
exchange integrals shall be restricted within a tight-binding 
approximation to nearest neighbour coupling, 



T. 



a/3 _ r«/3 



J i,j+A 



ja/3 ra/3 jaa . ra,/3±l jaf3 

ij ~ °i,j+A J "T °ij J ' 



(50) 
(51) 



respectively. Here A denotes the relative positions of nearest 
neighbours, both within the same layer, for s.c. (100): A = 
(0,1), (0,1), (1,0), (1,0). T al3 is the hopping between the 
layers a and (3 = a ± 1 and J aa is the exchange interaction 



within the layer a. For the following discussion, furthermore 
the hopping integrals and the exchange interaction have been 
assumed to be uniform within the whole film, 



rpaOi rj-iCt/3 rp 



= = Jff, 



(52) 



where Jff should not be mixed up with the s-f exchange in- 
teraction J. For explicit values we choose T = — O.leV, 
Jff = O.OleV, and the single-ion anisotropy Do/ Jf f = 0.01. 



A. /-spin correlation functions 

Before we can start calculating the electronic excitation 
spectra, we first have to evaluate the local-moment system 
considered in Sec. II B. 




FIG. 1. Layer-dependent magnetizations (S^) of a 20-layer 
s.c. (100) film for the layers a — 1, 2, . . . , 10. The magnetizations 
increase from the surface layer towards the center of the film. Inset: 
Dependence of the Curie temperature Tc on the film thickness n of 
s.c.-(lOO) films. 




T / T„ 



FIG. 2. The different /-spin correlation functions from Eqs. ( |33| ) 
calculated for the centre layers (a = 10, 11) of the 20-layer-film 
from Fig. |IJ 
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For example, Fig. [T] shows the layer-dependent magnetiza- 
tions of a 20-layer s.c.(lOO) film as a function of temperature. 
As for all other film thicknesses the layer-dependent magne- 
tizations increase from the surface layers (a = 1,20) 
towards the center layers {a = 10, 11) of the film. The inset 
of Fig. [j] displays dependence of the Curie temperature on the 
film thickness n. The complete set of /-spin correlation func- 
tions according to Eqs. (Q , (Q, (Q), and (f|9|) calculated for 
the center layer of the 20-layer film can be seen in Fig. ^[ 
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FIG. 3. Spectral density Ski(E) of a s.c.-(100)-monolayer for 
ferromagnetic saturation of the local-moment system (T = 0, 
(S z ) — S) as a function of energy and wavevector from T = (0, 0) 
to M = (tv, 7r) for T = and different values of the s-f exchange 
interaction J. 



B. The temperature-dependent electronic structure 

We discuss our results in terms of the spectral density 
S^"(E), defined in Eq. (|l2|), and the local quasiparticle den- 
sity of states, Eq. ( 13 ). We start our discussion of the temper- 



ature dependent electronic bandstructure with a special limit- 
ing case which gives us an insight into the underlying physics 
of the problem. The special limiting case of ferromagnetic 
saturation, T = 0, and empty conductionpband_-p = 0, is ex- 
actly solvah 



geometne: 



both, for the bulk materiaOE3El and for film 
The limiting case, therefore, provides a good 
testing grou nd for the theory for finite temperatures presented 
in Sec. [I A . 

It turns out that for T = the a =|-spectrum is rather 
simple, since a |-electron has no chance to exchange its spin 
with the ferromagnetically saturated localized /-spin system. 
The quasiparticle bandstructure is therefore identical to the 
free Bloch dispersion, only rigidly shifted by a constant en- 
ergy amount of — h JS, due to the first Ising-like term in the 
s/Hamiltonian (§). 

Fig. U shows the spin-| spectral density of a s.c.-(lOO)- 
monolayer for the special case of ferromagnetic saturation, 
T = 0, for different s-f interactions J. For J = the spectral 
density represents a (5-function located at the point of the free 
two-dimensional Bloch dispersion. For small s-f exchange 
coupling, J > 0, a slight deformation of the original Bloch 
dispersion sets in and the quasiparticle peaks get a finite width 
indicating a finite lifetime. For intermediate and strong cou- 
plings the spectral density splits into two parts corresponding 
to two different spin exchange processes between the excited 
spin-| electron and the localized /-spin system. The higher 
energetic part of the spectrum represents a polarization of the 
immediate spin neighbourhood of the electron due to a re- 
peated emission and reabsorption of magnons. The result is 
a polaron-like quasiparticle called the magnetic polaron. The 
low-energetic part of the spectrum is a scattering band which 
corresponds to the simple emission of a magnon by the spin- J, 
electroiy^yhich is necessarily connected with a spin-flip of the 
electrorcJ. 




J=0 
J=0.1 



J=0.2 



J=0.3 



FIG. 4. Density of states px(E) of a s.c.-(100)-monolayer for 
T = and for different s-f interactions J. 
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From the spectral density of Fig. || we get, using Eq. (|l3|), 
the local quasiparticle density of states pi (E) of a monolayer, 
displayed in Fig. |j. Here we see that the splitting of the spec- 
tral density discussed above transfers itself to the quasiparticle 
density of states as a gap for J 0.2. As for the spectral den- 
sity, the density of states of the spin-] electron is only rigidly 
shifted and therefore not displayed. 

However, this does not hold any longer for finite tempera- 
tures, T > 0. Fig. H exhibits the density of states of a s.c- 
(lOO)-monolayer for different s-f interactions and different 
temperatures. The dotted lines represent the case of vanish- 
ing s-f exchange, J = 0, where spin-] and spin-] spectra are 
equal. Since the electrons are not coupled to the local-moment 
system, we also have no temperature dependence. For finite 
s-f interaction we see from Fig. ^ that in the spin- J, density 
of states spectral weight is transferred from the high-energetic 
polaron peak to the low-energetic scattering peak. To explain 
this effect we have to consider the elementary processes which 
build the spectrum. The low-energetic scattering peak of the 
spin-], electron consists of two elementary processes. 



a=T J = 0,0.1,.. .,0.4 ^ 




-1 E / eV 1 



FIG. 5. Density of states p a (E) of a s.c.-(100)-monolayer for dif- 
ferent s -/interactions J and different temperatures T (in units of To). 
The dotted lines represent the case of J = where there is no dis- 
tinction between spin- J, and spin-] electron, Pi(E) = p^(E). 
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FIG. 6. Spectral density, S ka (E) = S&(E) = S&(E), of a 
s.c.-(lOO) double-layer as a function of energy and wavevector from 
r = (0, 0) to M = (-7T, 7r) for J = 0.2eV and different temperatures 
T (in units of Tb). 

Because of finite deviation of the /-spin system from sat- 
uration for T > 0, the .[-electron has a finite probability of 
entering the local frame as spin-] electron. This probability 
is zero for (S z ) = S (T = 0) and increases with increasing 
temperature. On the other hand, the spin- j electron can first 
emit a magnon and by that process reverse its spin, becoming 
a spin-f electron in the external frame of coordinates. The 
spectral weight produced by the first elementary process re- 
duces the spectral weight of the high-energetic polaron peak 
therefore shifting spectral weight from the high-energetic po- 
laron peaks towards the low-energetic scattering peak. 

For the spin-| electron we see from Fig. [| that for finite 
temperatures an additional peak rises at the high-energetic 
side of the spectra with increasing temperature. We can ex- 
plain this effect by the spin-| electron absorbing a magnon 
and subsequently, as spin-] electron forming a polaron. Here 
the magnon absorption by a spin-] electron is equivalent to 
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the magnon emission by a spin- j electron. In the case of ferro- 
magnetic saturation the system does not contain any magnons, 
which is the reason why there is no scattering peak in the spin- 
| spectrum at T = 0. As a result of the shifting of spectral 
weights towards lower energies for the spin-| electron and to- 
wards higher energies for the spin-| electron the densities of 
states for the two spin directions approach each other with in- 
creasing temperature. 

T/T c = 0, 0.4, 0.7, 0.9, 0.98, 1 
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J=0.2 
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FIG. 7. Density of states, p a {E) = pl(E) = pl(E), of a 
s.c.-(lOO) double-layer for J = 0.2 and different temperatures T 
(in units of Tc)- The curves for T — are furthest away from the 
paramagnetic solution, T = Tc (fat line). 



In the limiting case of T — > Tc the system has eventually 
lost its ability to distinct between the two possible spin direc- 
tions of the test electron because of the loss of magnetization 
of the underlying local-moment system, (S z ) — > 0. Hence as 
for the case of vanishing s-f interaction for T — Tc the den- 
sity of states of states of the spin-| electron equals that of the 
spin-| electron. Another feature which can be seen from Fig. 
| is that the positions of the four quasiparticle subbands, two 
for each spin direction, do not change with temperature. 

To further discuss the temperature effects we present with 
Figs. U and [7] the spectral density and the local density of 
states, respectively, of a s.c.-(lOO) double-layer (n = 2) for 
J = 0.2 and different temperatures. Again we see that the 
spectra for the two spin directions approach each other for 
T — > Tc- Another feature which can already be observed in 
Fig. U is that the increase of temperature results in the narrow- 
ing of the subbands. For the case of intermediate coupling, 
J = 0.2, according to Figs. || and ^ this band narrowing re- 
sults in the opening of a gap between the scattering and the 
polaron band with increasing temperature. 

This temperature enhanced band splitting has already been 
found for the three-dimensional caseEl It can be explained for 
the spin-| electron by the fact that for propagating in its own 
low-energetic subband it needs to find an appropriate lattice 
site. 
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FIG. 8. Local density of states, pf(E), of the first (a = 1, 20), 
second (a = 2, 19), and center (a = 10, 11) layer of a 20-layer 
s.c.-(lOO) film for different s-f interactions J and different temper- 
atures T/Tc = 0, 0.4, 0.7, 0.9, 0.98, 1. The curves for T = are 
furthest away from the paramagnetic solution, T = Tc (fat line). 

In the case of ferromagnetic saturation there is no restric- 
tion for the propagation of the spin-| electrons since spin-flip 
processes are impossible. With increasing temperature there 
is an increasing deviation of the local moments resulting in 
the possible magnon absorption by the spin-| electron and 
subsequent changing to the higher energetic polaron subband. 
Hence, the spin-| electron, to propagate in its own subband 
needs to move further distances with increasing temperature 
resulting in a reduced effective hopping and in a decreased 
bandwidth. 
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FIG. 9. Spectral density, , as a function of wavevector k 
from T = (0, 0) to M = (71", n) (horizontal axes) and energy E 
from -leV to leV (vertical axes) of the first (a = 1,20), second 
(a = 2, 19), and center (a = 10, 11) layer of a 20-layer s.c.-(lOO) 
film for J = 0.2eV and different temperatures (in units of Tc). 

In addition to the discussed temperature effects, Figs. |6| and 
[7] exhibit a typical two-peak structure which is caused by the 
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coupling of the two layers. This two-peak structure is replaced 
in the case of a film consisting of n equivalent layers by an n- 
peak structure. Generally, the spectra of the discussed local- 
moment films are characterised by an interplay between cor- 
relation (J), temperature effects, and geometry of the film. 
Figs. H and ^ display results for the local density of states and 
the layer-dependent spectral density of a 20-layer s.c.-(lOO)- 
film. Additionally to the dependence on the s-f exchange in- 
teraction and the temperature dependence we notice that the 
spectral density and the density of states show a typical layer- 
dependence due to the broken translational symmetry at the 
surfaces of the filmcii For the centre layers (a — 10,11) 
of the 20-layer- film we see from Fig. || that the local density 
of states of the spin-| electron at T = has already become 
pretty similar to the well-known tight-binding density of states 
of the three-dimensional s.c. lattice whereas the density of 
states of the surface layers (a = 1, 20) exhibits the character- 
istic semi-elliptic profile. 

IV. SUMMARY 

We have investigated the electronic quasiparticle bandstruc- 
ture of a ferromagnetic semiconductor film. A single test elec- 
tron (s-band) is coupled by an intra-atomic s-f inter-band ex- 
change to a system of localized 4/-moments. That may be 
regarded as a proper model description for Euo and EuS. Our 
approach uses a moment-conserving decoupling procedure for 
suitable defined Green functions. The fact that our theory 
evolves continuously from the exactly solvable limiting case 
of ferromagnetic saturationES gives it a certain trustworthi- 
ness. 

The exchange coupling of the conduction electron to the 
local-moment system gives rise to a correlation induced split- 
ting of the quasiparticle spectra. A polaron part may be in- 
terpreted as a repeated emission and reabsorption of magnons 
by the conduction electron resulting in a new quasiparticle, 
the magnetic polaron. A rather broad scattering peak is due 
to a simple magnon emission or absorption by the conduc- 
tion electron. This pronounced splitting depends on the ac- 
tual value of the exchange interaction J. For small values of 
J only a renormalization of the one-electron energy occurs 
resulting in a deformation of the free Bloch dispersion. For 
higher values of J, the mentioned splitting of the spectra into 
polaron part and scattering part sets in. 

We intend to apply the presented model to study the tem- 
perature dependent electronic structure of EuO and EuS films. 
Therefore the electronic part of the Hamiltonian (^J) and the s- 
f exchange interaction (H) have to be modified to include the 
multiband aspect of real substances. This will be done by sub- 
stituting the tight-binding bandstructure ( p0[ ) by a realistic one 
taken from a bandstructure calculation. 

Another highly interesting field which we want to use our 
theory for is the evaluation of the temperature dependence of 
surface states. In a previous paper we have calculated surface 
states for the special case of T = and n-^ by modify- 
ing the hopping in the vicinity of the surfaceES. The extension 



of these calculations to finite temperatures promises to give 
an understanding of recent experimental results concerning 
the temperature dependence of electronic states on surfaces 
of rare earths. 
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